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ABSTRACT 

The starburst galaxy NGC 253 was observed with the Energetic Gamma Ray Experiment Telescope 
(EGRET) aboard the Compton Gamma Ray Observatory (CGRO) satellite. We obtain a 2 o upper limit to 
the y-ray emission above 100 MeV of 8 x 10 8 photons cm 2 s . Because of their large gas column 
densities and supernova rates, nearby starburst galaxies were predicted to have y-ray fluxes detectable by 
EGRET. Our nondetection of y-rays from NGC 253 motivates us to reexamine in detail the premise of 
supernova acceleration of cosmic rays and the effect of enhanced cloud densities, photon densities, and 
magnetic fields on the high-energy spectra of galaxies. By modeling the expected y-ray and synchrotron 
spectra from NGC 253, we find that up to 20% of the energy from supernovae is transferred to cosmic 
rays in the starburst, which is consistent with supernova acceleration models. Our calculations match the 
EGRET and radio data well with a supernova rate of 0.08 yr *, a magnetic field B > 5 x 10 5 G, a 
density n ~ 300 cm \ a photon density U ph ~ 200 eV cm 3 , and an escape timescale r 0 < 10 Myr. 
Subject headings: cosmic rays galaxies: individual (NGC 253) — galaxies: starburst 
gamma rays: observations gamma rays: theory stars: statistics 


1. INTRODUCTION 

1.1. Diffuse Gamma-Ray Emission in Spiral Galaxies 

Since the earliest high-energy experiments, the diffuse 
y-ray emission in the Milky Way has been shown to corre- 
late closely with gas column density (see review by Bloemen 
1989). In the interstellar medium (ISM), high-energy y-rays 
(E > 100 MeV) are produced mostly through interactions 
between matter and cosmic-ray (CR) nuclei (inelastic 
collisions), and from electrons (bremsstrahlung emission 
and inverse Compton scattering). Cosmic rays are thought 
to originate in the shock fronts of supernova explosions 
(e.g., Ginzburg & Syrovatskii 1964). Therefore, the y-ray 
flux is a measure of the supernova rate and the average gas 
column density encountered by a CR before it escapes a 
galaxy. Hence, galaxies with high supernova rates and large 
gas masses should have more high-energy emission. 

Of all galaxies, starburst galaxies have the highest 
inferred supernova rates and are predicted to have CR den- 
sities several hundred times larger than the Milky Way 
(Suchkov, Allen, & Heckman 1993). The large far-infrared 
(FIR) luminosity of a starburst galaxy implies that the rate 
of massive star formation is very high, since this emission is 
presumed to be starlight reradiated by dust (Rieke et al. 
1980). Starbursts also contain large masses of molecular gas 
as traced by CO (Young et al. 1989), and a significant frac- 
tion of that gas is dense, with > 10 4 cm 3 (Paglione et 
al. 1995; Jackson et al. 1995; Heifer & Blitz 1993; Wild et al. 
1992; Sage, Shore, & Solomon 1990; Carlstrom 1989; 
Mauersberger & Henkel 1989). 

1.2. The Starburst Galaxy NGC 253 

NGC 253 is one of the archetypal starburst galaxies iden- 
tified by Rieke et al. (1980). A supernova rate of roughly 0.1 
has been derived from high-resolution radio contin- 
which show many nonthermal point sources 
^■'''AJlyestad 1988). From CO observations, the 


gas column density of the inner kiloparsec of NGC 253 is 
estimated to be ~10 23 cm 2 (Canzian, Mundy, & Scoville 
1988), roughly 10 times that of the Milky Way (e g., Sanders, 
Solomon, & Scoville 1984). From multiwavelength molecu- 
lar line studies, the average density of the molecular clouds 
in NGC 253 is roughly 10 5 cm -3 , with a volume filling 
factor of < 10 3 (Paglione et al. 1995; Jackson et al. 1995; 
Wall et al. 1991; Mauersberger & Henkel 1989). The 
average gas density over the starburst volume is therefore 
30 300 cm 3 , up to 100 times that of the Milky Way. 
(Throughout this paper we assume that the distance to 
NGC 253 is 3.4 Mpc.) 

The continuum spectrum of NGC 253 peaks in the FIR 
at ~ 100 /xm with a luminosity of L F , R = 3 x 10 10 L Q (Rice 
et al. 1988), roughly 5 times the FIR luminosity of the Milky 
Way (Cox & Mezger 1989). Because the IR region is con- 
centrated within the central 500 pc of the starburst nucleus 
(Pina et al. 1992; Heckman, Armus, & Miley 1990), the 
photon density in this region of NGC 253 is much larger 
than in the Milky Way. Therefore, inverse Compton scat- 
tering of the FIR radiation off the large population of CR 
electrons in NGC 253 (Carilli et al. 1992) should be signifi- 
cant (Goldshmidt & Rephaeli 1995). 

Another well-known starburst galaxy, M82, which has 
similar gas properties and nearly the same supernova rate 
as NGC 253, may have a CR density roughly 300 times 
higher than that of the Milky Way (Suchkov et al. 1993). 
Because the y-ray flux from M82 was predicted to be detect- 
able with EGRET (Akyiiz, Brouillet, & Ozel 1991), we 
anticipated that NGC 253 would be detectable in hard 
y-rays as well. However, preliminary EGRET results show 
only upper limits to the y-ray emission from M82 and NGC 
253 (Sreekumar et al. 1994). Here we analyze EGRET 
observations of y-ray emission from NGC 253 and examine 
in detail the high-energy flux expected from a starburst 
galaxy. 
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2. OBSERVATIONS AND RESULTS 

The EGRET instrument is sensitive to photons in the 
energy range from roughly 20 MeV to 30 GeV. A complete 
description of EGRET is given by Hughes et al. (1980) and 
Kanbach et al. (1988, 1989). Thompson et al. (1993) present 
the results of the instrument calibration. 

NGC 253 was observed for a total on-axis integration 
time of roughly 2 weeks. We subtracted a background 
modeled after the Galactic H i emission. This measurement 
is very sensitive, since NGC 253 is located at the south 
Galactic pole where the background is minimal. We thereby 
obtain a 2 a upper limit of 8 x 10 8 photons cm 2 s 1 for 
the flux above 100 MeV from NGC 253. This limit we find 
from reanalyzing the EGRET data is slightly different from, 
yet consistent with, that found by Sreekumar et al. (1994). 

3. THE COSMIC-RAY ELECTRON DISTRIBUTION 
IN CLOUDS 

The lack of significant y-ray emission from M82 is consis- 
tent with the predictions of Akyiiz et al. (1991) and Pohl 
(1994), given the sensitivity of EGRET. However, this result 
does not immediately explain our nondetection from NGC 
253 since the smaller size, higher average molecular gas 
density, and younger age of the starburst in NGC 253 imply 
that the CR density (n CR ) and y-ray emissivity of NGC 253 
should be higher than that of M82. The starburst volume of 
NGC 253 is smaller than that of M82 in radio, infrared, and 
CO maps (Pina et al. 1992; Telesco et al. 1991 ; Carlstrom et 
al. 1989; Canzian et al. 1988), which indicates not only an 
increase in the CR density in NGC 253 despite a lower 
supernova rate, but also an increased photon density, which 
yields a higher inverse Compton flux. The larger average 
molecular gas density in NGC 253 (Jackson et al. 1995) 
implies higher secondary electron and neutral pion pro- 
duction, which in turn lead to enhanced y-ray emission. 
Rieke, Lebofski, & Walker (1988) postulated that the star- 
burst in NGC 253 is less evolved than that of M82. Also, 
radio spectral index maps imply that diffusion rather than 
convection is still a major escape mechanism for CR elec- 
trons in NGC 253 (Carilli et al. 1992). Therefore, cosmic 
rays should remain near molecular clouds longer in NGC 
253 than in M82 (see Seaquist & Odegard 1991) and 
thereby support reactions that lead to y-ray emission. 
Finally, several other independent methods for estimating 
n CR in starburst galaxies suggest that the CR density is 
higher in NGC 253 than in M82 by factors of a few 
(Suchkov et al. 1993). 

Therefore, we proceed to calculate the expected y-ray flux 
from NGC 253 based on the observed supernova rate, gas 
density, and synchrotron radio emission. 

3.1. Modeling the Interstellar Cosmic-Ray Distribution 

To calculate the CR density and y-ray flux from a star- 
burst, we model the CR production and radiation mecha- 
nisms from a region with dense clouds and a high 
supernova rate. An initial CR spectrum must be assumed, 
and the energy distribution of CR electrons must be gener- 
ated. We consider the contribution of both primary and 
secondary electrons, since secondary production in the 
dense environment of a starburst nucleus is probably sig- 
nificant (Marscher & Brown 1978). We also examine the 
effects of a large photon density, magnetic field, and cloud 


density on the CR energy distribution and subsequent y-ray 
emissivity. 

In § 3.2 we outline the procedure for generating the 
steady state CR energy distribution in the ISM. The dis- 
cussion concentrates on secondary electron production and 
loss mechanisms. We discuss the final form of the secondary 
electron distribution in § 3.3, and we describe the primary 
CR electron spectrum in § 3.5. We determine the input CR 
proton spectrum in § 3.4. We discuss the consequent radio 
and y-ray emission in § 4. 

3.2. Solution to the Diffusion-Loss Equation 

The procedure we follow for generating the secondary 
electron spectrum of an interstellar cloud was outlined by 
Marscher & Brown (1978). The general steady state 
diffusion-loss equation for CR electrons may be written as 
(e.g., Longair 1981) 

N(E) d , 

^ — [b(E)N(E)-] - Q(E) = 0 . (1) 

The electron energy distribution (the number of electrons 
with energy between E and E + dE) is N(E). The source 
function Q(E ) is the sum of the production rates of second- 
ary electrons. The total time-independent loss rate is 1 /t(E). 
The electron energy loss rate, b(E) = — dE/dt , is given by 
(e g., Ginzburg & Syrovatskii 1964) 

b(E) = 2.52 x 10 17 n[6.27 + ln(£/mc 2 )] 4- (V/R)E 

+ E02 x 10 16 [U ph + 2.49 x 10 lo £ 2 ]£ 2 GeV s 1 . 

( 2 ) 

Ionization losses, which are important at lower energies, are 
represented in the first term of equation (2), where n is the 
average interstellar density, £ is the particle energy, and me 2 
is the electron rest energy. The second term is the loss rate 
due to adiabatic expansion. For NGC 253, the starburst 
wind speed is V ~ 300 km s 1 and the starburst radius is 
R = 325 pc (Heckman et al. 1990; Canzian et al. 1988; 
McCarthy, Heckman, & van Breugel 1987). Adiabatic 
expansion losses have a minor effect on the electron dis- 
tribution and may be treated as an effective sink. The last 
term is the loss rate due to synchrotron radiation and 
inverse Compton scattering, where the photon energy 
density U ph is in eV cm -3 , and the magnetic field B is in 
gauss. 

We assume that the change in the electron distribution 
due to bremsstrahlung radiation may be treated with a 
single characteristic loss rate 1 /t b , and that an electron loses 
all its energy in the interaction. The total time-independent 
loss rate is comprised of the energy-dependent diffusion rate 
and the bremsstrahlung loss rate (Blandford & Ostriker 
1980; Ormes & Freier 1978), 

r *( E) = x b \E) + t d \E ), 

where 

r B l (E) = 9.6 x 10 l6 n s 1 (3) 

and 

T D '(E) = T 0 'fiE 0 Ge \ + T c . 

The characteristic escape time in a “ leaky box ” model is t d , 
where is the particle velocity in light units, and x c is the 
convection timescale. We assume that the convection-loss 
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timescale is within an order of magnitude of the energy- 
dependent escape time t 0 at 1 GeV based on the evidence of 
supernova winds in the nucleus of NGC 253 (McCarthy et 
al. 1987; Antonucci & Ulvestad 1988) and convective 
“chimneys” in its disk (Carilli et al. 1992). Our calculations 
include the cases t 0 = 0.1, 1, and 10 Myr, and we assume 
t c « 1 Myr. 

Equation (1) may be solved analytically with the inte- 
grating factor 

W£)e,p [“£5^5i]' (4) 

The equation for the electron energy spectrum in a cloud is 
then 


N(E) = Q(E') 


P Q{E’) 

Jfo 


dE 


dE 


exp [ — F(E\ E)~\dt f , 


( 5 ) 


where F(E\ E) is defined as the integral in equation (4). 


3.3. Secondary Electron Production 

The three reactions between CRs and interstellar matter 
that produce most of the interstellar secondary electrons are 
knock-on collisions, neutron beta decays, and charged-pion 
decays. Each source term can be written as 


Qi(E) oc n 


m*- 


( 6 ) 


where dJ p /dE p = N p (E p )Pc/4n is the differential energy spec- 
trum of CR protons, and da(E , E p )/dE p is the differential 
cross section for the interaction. The CR proton spectrum 
N p will be discussed in § 3.4. 

Knock-on or Coulomb collisions produce CR electrons 
when a high-energy particle transfers energy much greater 
than the binding energy to an atomic electron. We use the 
form of the knock-on cross section given by Abraham, 
Brunstein, & Cline (1966). The resulting spectrum of sec- 
ondary electrons follows an approximately constant power 
law that is up to 20% steeper than the original proton 


injection spectrum. The secondary electron spectrum from 
neutron beta decays is roughly an order of magnitude lower 
than that from knock-on collisions and has a similar slope 
(Abraham et al. 1966). 

The above two mechanisms produce most of the 
lower energy CR secondaries. Charged-pion decay 
(tt* -» ► /i* -> e*) dominates the production of higher energy 
CR electrons (E > 100 MeV). Dermer (1986) provides a 
concise overview and comparison of past studies of pion 
production in the ISM. We follow the procedure outlined 
by Marscher & Brown (1978). The pion and subsequent 
electron source spectra peak below 100 MeV (&^m K c 2 ), 
and the slope approaches that of the incident CR proton 
spectrum above 1 GeV (see Dermer 1986). 

The total secondary electron distribution is the summa- 
tion of the above source spectra. The result is displayed in 
Figure 1. The steady state secondary CR electron distribu- 
tion calculated from the total source spectrum with equa- 
tion (5) is shown in Figure 2. The high-energy electron 
spectrum is steeper for an increased magnetic field and 
photon field owing to enhanced synchrotron and inverse 
Compton losses. Bremsstrahlung and adiabatic expansion 
losses decrease the entire spectrum while maintaining the 
shape and are most pronounced, yet still a minor effect, near 
1 GeV. A shorter characteristic escape time decreases the 
high-energy CR spectrum but does not change the slope. 
Escape losses have a stronger effect on the total production 
rate of secondary electrons, since these losses have a strong- 
er bearing on the steady state CR proton distribution (see 
§ 3.4). 


3.4. Cosmic-Ray Proton Energy Spectrum 

The proton energy distribution is required to derive the 
secondary electron source functions in § 3.3. The proton 
spectrum is usually assumed to follow a power law pro- 
portional to E p s , where s = 2.5 2.75 approximates the local 
CR distribution (Dermer 1986). However, because the goal 
of this project is to test the mechanism of supernova acceler- 
ation of cosmic rays, we constrain only the CR injection 
spectrum, and we assume that it is constant with time. The 



Fig. 1. The total secondary CR electron source spectrum ( heavy line) from knock-on (Q t ), neutron beta decay and pion-muon-electron decay {Q K(ie ) 
reactions ( dashed lines) for two values of the density n and escape time x 0 in a “ leaky box ” model of diffusion. 
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Fig. 2. The steady state primary CR electron spectrum ( continuous line) compared with that of secondary electrons ( dashed line) for n = 100 cm 3 
B = 5 x 10 5 G, U ph = 200 eV cm 3 , t 0 = 10 Myr, and N p /N e = 100. 


initial proton spectrum N p (E p ) takes the form of equation 
(5) with an injection spectrum Q cc E p \ This CR distribu- 
tion changes with time owing to ionization and escape 
losses as the cosmic rays encounter the ISM. The loss time 
t (E p ) is now given only by t d in equation (3), and the energy 
loss rate is (e.g., Ginzburg & Syrovatskii 1964) 


b(E p ) — 1.84 x 10 


7 % ll.l+21ny, + ln0 



-I- 0 .65nca n E p GeV s 1 , (7) 

where v is the particle velocity, and y p is the total proton 
energy divided by the proton rest energy Me 2 . For protons 
with total energy > 1.22 GeV, losses due to pion production 
must also be included. For this term, we use the relation 


given by Mannheim & Schlickeiser (1994), where a K = 3 
x 10" 26 cm 2 is the total pion cross section, and E 
+ Me 2 > 1.22 GeV. P 

The value of the spectral index s is estimated from pre- 
vious observations of Galactic and extragalactic super- 
novae (Reynolds & Ellison 1992; Huang et al. 1994), as well 
as the local Galactic CR distribution (e.g., Ormes & Freier 
1978) and the nucleus of M82 (Pohl 1994). Indices between 
2 and 2.6 are observed for most supernova remnants and 
are predicted from blast-wave theory and models of first- 
order Fermi acceleration of CRs (Bell 1978a, b; Reynolds & 
Ellison 1992, and references within). We consider two 
models with s = 2.2 and 2.5. 

The calculated steady state proton spectrum is shown in 
Figure 3. The distribution with s = 2.2 compares well with 



Fig. 3. - The steady state CR proton spectrum for two values of the density and escape timescale. The kink just below 0.3 GeV results from the threshold 
to pion production. 
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that of the Milky Way (see Dermer 1986; Ormes & Freier 
1978; Bell 1978a, b). Note that the spectrum decreases with 
density at lower energies because of ionization losses. The 
high-energy protons reach a steady state by either the char- 
acteristic escape time t 0 , or the time for energy loss from 
pion production t„ ~ 0.5(n/100 cm 3 ) 1 Myr. Therefore at 
higher densities, pion production dominates the energy 
losses above 1 GeV, but for lower densities, the pion pro- 
duction and escape loss rates are comparable (Fig. 3). 


3.5. Primary Electron Distribution 

The contribution of primary CR electrons to the diffuse 
y-ray emission in NGC 253 must also be determined. 
Primary electrons are presumably accelerated in supernova 
explosions in the same way as protons; the injection spec- 
trum is a power law with s = 2-2.6. Figure 2 shows the 
steady state primary electron distribution calculated from 
equation (5). 


3.6. N ormalization of the CR Injection Spectrum 

In § 3.4 we introduced the CR proton injection spectrum 
as Q = KE s . The constant of proportionality K may be 
obtained from the total power transferred by supernovae 
into cosmic rays within a given volume 



QEdE = riP'V/r . 


( 8 ) 


Here rjP is the energy transferred from a supernova to CRs, 
where P « 10 51 ergs, and r} is the transfer efficiency. The 
supernova rate is *P, and tC is the volume. Solving for the 
normalization, we find (for s > 2) 


riP'V , 

k = - l —(s 
i- 


2 )E’ mla 2 . 


(9) 


We treat the case of s = 2.2 and E min = 100 MeV. 

An additional scaling factor, the proton-to-electron flux 
ratio N p /N e , must be included for primary electrons. This 
ratio is 100 for the Milky Way, but it could be smaller if 
CRs are accelerated by multiple shocks (Bell 1978b). X-ray 
data for M82 are consistent with N p /N„ as low as 30, 
although a ratio of 100, which implies a lower energy trans- 
fer efficiency, may be more appropriate (Volk, Klien, & 
Wielebinski 1989a). 


4. DISCUSSION 

We use the calculations of § 3 to model the radio and 
y-ray emission from the starburst galaxy NGC 253. The 
starburst nucleus is assumed to be a uniform cylinder with a 
constant supernova rate. We have examined the expected 
emission from the disk of NGC 253 as well, where the 
source of CRs is the population leaving the nuclear region. 
However, we find the high-energy emission from the disk to 
be negligible, and the parameters needed for the calcu- 
lations are poorly constrained by observations. Therefore, 
we approximate that all the y-ray emission from NGC 253 
originates in the nucleus, and once a cosmic ray leaves the 
starburst region, it is no longer considered in the model. 

4.1. Radio Emission 

Radio images of NGC 253 show a prominent nucleus 
which matches the dimensions of the CO and IR starburst 
regions (Pina et al. 1992; Antonucci & Ulvestad 1988; 


Canzian et al. 1988). The spectrum from the nucleus is flat 
(Ulvestad & Antonucci 1991), and there is also considerable 
nonthermal emission within the disk and in a large, steeper 
spectrum, diffuse halo around the galaxy (Beck et al. 1994; 
Carilli et al. 1992; Klein et al. 1983; Beck et al. 1979). 

We calculate the synchrotron emissivity e v with the elec- 
tron distribution from § 3 assuming an isotropic distribu- 
tion of pitch angles. We include thermal free-free emission 
// with T = 8000 K, as well as free-free absorption and dust 
emission in order to determine the total continuum flux 
density of NGC 253. We vary the emission measure EM so 
that the thermal flux density matches millimeter obser- 
vations (Carlstrom et al. 1989; Elias et al. 1978). The total 
radio continuum flux density is then given by 

S r =fJn + ^-l\ -e " u )Ql , (10) 

T v 

where the free-free optical depth i" oc EM. The source solid 
angle ft - 300 arcsec 2 and the line-of-sight depth / - 630 pc 
are estimated from CO interferometry and 6 cm images of 
NGC 253 (Canzian et al. 1988; Antonucci & Ulvestad 
1988). We also calculate the dust emission assuming 
S v (dust)oc v i 5 B v (T D \ where B v is the Planck function and 
T d * 40 K (Telesco & Harper 1980). 

For normalization, we choose a cylindrical starburst 
region 70 pc thick, with a radius of 325 pc (Canzian et al. 
1988; Antonucci & Ulvestad 1988), and T* = 0.08 yr 1 (Van 
Buren & Greenhouse 1994; Ulvestad & Antonucci 1994). 
The total modeled radio and FIR spectrum from the 
nucleus of NGC 253 is shown in Figure 4. The model repro- 
duces the radio spectrum of the nucleus of NGC 253 for the 
following values of the free parameters: s = 2.2, B>5 
x 10 5 G, n = 300 cm -3 , U ph = 200 eV cm 3 , t 0 = 10 
Myr, and EM ~ 10 5 cm 6 pc. 

If we consider only the emission from secondary elec- 
trons, the required efficiency of energy transfer from super- 
novae to CRs is roughly 50%. Either CR acceleration in 
starburst galaxies is very efficient, or primary electrons con- 
tribute in part to the emission. To include emission from 
primary electrons, we must evaluate the proton-to-electron 
ratio, and the choice of this ratio affects the efficiency. 
Figure 5 shows how this ratio N p /N e depends on the energy 
transfer efficiency, t]. Also shown is the fraction of the total 
flux due to emission from secondaries F 2dy /F Xol . Given 
N p /N e = 100, we find rj ~ 20%, which is somewhat higher 
than values predicted for the Milky Way (see Volk et al. 
1989a; Volk, Klien, & Wielebinski 1989b; and references 
therein). In the case of a lower N p /N e ratio of 30, which is 
possible if CRs are shock accelerated more than once (Bell 
1978b), the efficiency necessary to match the radio emission 
is only 9%. Note that the contribution from secondary elec- 
trons is significant in a starburst nucleus, whereas primary 
electrons dominate Galactic synchrotron emission 
(Marscher & Brown 1978). 

4.2. Gamma-Ray Emission 

4.2.1. Cosmic-Ray Electrons 

Bremsstrahlung emission and inverse Compton scat- 
tering of CR electrons are prominent sources of soft y-rays 
(E < 100 MeV). In the Galaxy, inverse Compton scattering 
is not so important (Bloemen 1989), but it may be in star- 
burst galaxies because of their high photon densities. Simi- 
larly, the bremsstrahlung emission may be enhanced owing 
to the higher average densities. 
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F>°- 4. — The radio and FIR spectrum of NGC 253. The triangles are the integrated flux densities (Carlstrom et al. 1989; Klein et al. 1983; Ruhr et al. 
1981 ; Telesco & Harper 1980; Elias et al. 1978). The other symbols are the intensities from only the nucleus ( crosses , Carilli et al. 1992; squares, Hummel, 
Smith, & van der Hulst 1984; star, Klein & Emerson 1981 ; dots, Beck et al. 1979). The model synchrotron spectrum is the summation of the primary and 
secondary electron emission. 


We generate the inverse Compton emissivity with the CR 
electron distribution calculated in § 3 using well-known 
equations (Rybicki & Lightman 1979, p. 207). To model the 
emission from NGC 253, we again assume a photon energy 
density of U ph = 200 eV cm~ 3 from the FIR luminosity 
(Rice et al. 1988) and the size of the starburst region 
(Antonucci & Ulvestad 1988; Canzian et al. 1988; Pina et 
al. 1992). The bremsstrahlung emissivity is calculated from 
the CR electron spectrum with 

Q b (EJ = n<rcE y 1 N(E)dE cm" 3 s _1 GeV 1 . (11) 

jEy 

In the Galaxy, the contribution of inverse Compton 
y-rays to the total diffuse y-ray emission is roughly 10% of 
that from bremsstrahlung (e.g., Bloemen 1989). In NGC 
253, this ratio is preserved because the average cloud 


density and photon field (and consequently, the bremsstrah- 
lung and inverse Compton emission) are higher than in the 
Milky Way by roughly the same factor of ~ 100. 

4 . 2 . 2 . Neutral Pion Decays 

The dominant mechanism for generating diffuse y-ray 
emission (E > 100 MeV) in the ISM is neutral pion decays. 
The neutral pion spectrum is half the total charged pion 
spectrum discussed in § 3.3. Each n° decays into two y-rays. 
The source function of y-rays from n° decays is (Cavallo & 
Gould 1971) 


QAEy) = 2 f * 

jEi 


(i/gKygj 

(El ~ 4£g) 1 ' 2 


dE n cm 3 


s" 1 GeV 1 , 


where E l = E y + E 2 /E yy and E 0 = M n c 2 / 2. 


( 12 ) 
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Fig. 5. The proton to electron ratio NJN e and the fraction of emission from secondary electrons F 2dy /F toV vs. the efficiency rj by which the total energy 
of a supernova is converted into cosmic rays. 
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4.3. The Expected Gamma-Ray Spectrum 
from a Spiral Galaxy 

4.3.1. NGC 253 

The model y-ray spectrum of NGC 253 is shown in 
Figure 6. Above 100 MeV, the emission is dominated by 
neutral pion decays, and the flux near 1 MeV is due mostly 
to inverse Compton scattering. Our model greatly under- 
estimates the lower energy (0.05 MeV < E < 10 MeV) y-ray 
flux seen from OSSE observations (Bhattacharya et al. 
1994), but recent analysis of these data has shown that most 
of the emission arises not in the nucleus, but from the disk 
of NGC 253 (Goldshmidt & Rephaeli 1995). Our inverse 
Compton flux is somewhat lower than that modeled for the 
nucleus by Goldshmidt & Rephaeli (1995), but our results 
are not inconsistent with theirs considering the differences 
in our initial assumptions. The total integrated flux of y-rays 
above 100 MeV from Figure 6 is 4 x 10' 8 photons cm' 2 
s'. Moderate changes in the free parameters can decrease 
the model y-ray flux by up to a factor of 2 and still match 
the synchrotron spectrum. If we force our inverse Compton 
spectrum to match the OSSE data, the predicted flux above 
100 MeV is 5 x 10 8 photons cm -2 s \ which is still con- 
sistent with the EGRET limit. 

Typical age estimates for the starburst in NGC 253 
are ~ 10 7 yr (Rieke et al. 1980, 1988). The CR distribution 
should reach a steady state in this time. However, if it has 
not, the evolving CR spectrum is slightly different than the 
ones shown in Figure 3. Because the CR protons below 
1 GeV reach a steady state much more quickly than those at 
high energies, the evolving distribution resembles one in 
steady state spectrum with a lower escape timescale (see Fig. 
3). The resulting secondary electron and synchrotron 
spectra are slightly lower and flatter, and the n° decay y-ray 
flux is lower as well. The radio spectrum of NGC 253 is 
matched best by a steady state CR distribution given t 0 = 
10 Myr, and we conclude that the CRs have reached equi- 
librium. 

4.3.2. The Galaxy 

In comparison, the y-ray flux from the Galaxy (Fichtel et 
al. 1993; Bloemen 1989) is matched well by the model if we 


assume a supernova rate of 1/50 yr 1 (Tammann 1974) and 
a cylindrical volume 10 kpc in radius and 250 pc thick. The 
transfer efficiency necessary to match the COS B and 
EGRET data is t] ~ 5%, in good agreement with theoretical 
calculations (Markiewicz, Drury, & Volk 1990; Drury, 
Markiewicz, & Volk 1989). We assume B = 10 6 G, n = 1 
cm 3 , x 0 = 10 Myr, and V P * = 1 eV cm 3 in the Galaxy. 
Unlike the starburst galaxy, emission from primary elec- 
trons dominates the flux from the Milky Way, given 
N p /N e = 100. 

4.3.3. M82 

The synchrotron spectrum of M82 is flatter than that of 
NGC 253 (Klein, Wielebinski, & Morsi 1988), and it is 
matched by the model if we adopt somewhat lower values 
for the magnetic field, photon field, and characteristic 
escape time than those used for NGC 253. These changes 
are supported by the arguments at the beginning of § 3, and 
they result in a lower y-ray flux from M82, consistent with 
the EGRET limit of 5 x 10 8 photons cm' 2 s' 1 
(Sreekumar et al. 1994) and other model calculations (Pohl 
1994; Akyiiz et al. 1991). We should mention that the flatter 
synchrotron spectrum is also consistent with higher den- 
sities owing to the decrease in low-energy CRs through 
ionization losses. However, molecular line studies of these 
galaxies indicate that the average densities of the cloud 
cores in M82 are up to 10 times lower than those in NGC 
253 (Jackson et al. 1995), and higher densities also imply an 
increase in y-ray flux. Therefore, the average gas densities in 
M82 are most likely lower than those in NGC 253, and the 
flatter synchrotron spectrum of M82 is a result of a lower 
magnetic field and photon density. We cannot exclude the 
possibility that the CR production rates are different in 
these two galaxies either, since the physical conditions of 
the ISM in galaxies regulate both the production of y-rays 
and CR losses. Thus, changing the density and photon field 
tends to have little effect on the high-energy emission (Pohl 
1994). 

5. CONCLUSIONS 

We obtain a sensitive upper limit to the y-ray flux from 
the starburst galaxy NGC 253 of 8 x 10 8 photons cm 



Fig. 6. The model y-ray spectrum of NGC 253 from bremsstrahlung, inverse Compton, and n° decays 



302 


PAGLIONE ET AL. 


s" 1 with EGRET. For comparison, we determine the steady 
state energy distribution of CRs in the ISM of a starburst 
galaxy. The model fits the synchrotron spectrum from NGC 
253 well, given the high magnetic fields, gas densities, and 
photon fields in starburst nuclei. We derive a y-ray flux from 
the nucleus of NGC 253 of 4 x 10 -8 photons cm -2 s 1 . 
This result assumes a supernova rate of 0.08 yr -1 , a mag- 
netic field B > 5 x 10 5 G, an escape timescale t 0 ~ 10 
Myr, a density n ~ 300 cm" 3 , and a photon density U ph ~ 
200 eV cm 3 . From these calculations, we find that up to 
20% of the energy from supernovae in the starburst is trans- 
ferred to CRs, which is somewhat larger that the value 


derived for the Galaxy (3%— 10%), unless we assume a lower 
proton-to-electron ratio in NGC 253. 
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